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Abstract 

We report an investigation of laser frequency stabilization using a whispering gallery mode 
resonator that is temperature stabilized by a dual-mode technique. This dual-mode technique 
has yielded mode volume temperature instabilities at the nK level, suggesting that high fre- 
quency stability may also be reached. Here, we experimentally and theoretically investigate 
the dynamics of such a system and the important factors affecting the achievable frequency 
stability. We calculate that the dual-mode technique can reduce the effective fractional tem- 
perature coefficient of the reference system to 3.6 x 10~ 8 K _1 within the temperature feedback 
bandwidth. We demonstrate a 1560 nm laser stabilized to 1.3x 10~ 12 at 1 s and 1.1 x 10~ 10 at 
(— i 1 1000 s, corresponding to a long-term drift of 21 kHz/hr. 

1 Introduction 

Whispering gallery mode resonators (WGMRs) demonstrate extremely high quality factor and 
correspondingly narrow linewidth optical resonances over a large wavelength range, making 
them attractive for laser frequency stabilization applications. In fluoride resonators, qual- 
ity factor (0 in excess of one billion is regularly demonstrated, and Q >6xl0 10 has been 
achieved 0], corresponding to a linewidth of £=4.5 kHz - comparable to that of high-finesse 
Fabry Perot (FP) cavities. With special annealing techniques, an even higher Q of 3 x 10 ! 1 has 
been demonstrated (2). The wide transparency window of a crystalline WGMR mitigates the 
need for expensive mirror-coating runs and facilitates use of a single cavity at multiple wave- 
lengths, while their compact size and robustness to mechanical vibrations makes them well 
suited for applications requiring portable high-performance lasers. 

Impressive gains to laser linewidths and short-term stability have already been achieved using 
WGMRs. Linewidths of 13 kHz were demonstrated with a resonator acting as a transmission 
filter in a fiber laser cavity |3), and a 200 Hz linewidth was achieved through injection locking to 
a WGMR |4|. When employed as an external frequency reference cavity, crystalline resonators 
have delivered Allan deviations of 6x 10~ 14 at 100 ms 0. 

Frequency stability of WGMR reference cavities is limited by the fact that, in contrast to an 
FP resonator, the light travels in a solid medium. At short time scales, thermorefractive noise 
defines the limit and has been analyzed theoretically [6,7| and measured experimentally (8). 
At longer time scales, stability is limited by temperature changes in and around the resonator. 
These changes pull the frequency both through thermorefraction and thermal expansion, though 



the latter dominates in the common resonator host crystals magnesium fluoride and calcium 
fluoride. Ultra-low expansion (ULE) glasses have thermal expansion coefficients about three 
orders of magnitude lower than common optical crystals, but their poor transmission qualities 
preclude their use as a WGMR material. 

Recently, dual-mode temperature stabilization schemes have been investigated as a means to 
address this heavy dependence of frequency upon temperature l9 ffl"2l . In such schemes, the fre- 
quency spacing between modes with different temperature coefficients (tempcos) of frequency 
provides an extremely sensitive measurement of the mode volume temperature and facilitates 
active stabilization thereof. Mode volume temperature instabilities less than 10 nK have been 
demonstrated |9][10). But the WGMR represents a complex thermal system with several dif- 
ferent sources, sinks and the resulting gradients. We find that a highly stabilized mode volume 
temperature does not guarantee frequency stability, a finding also suggested by other recent 
work ifTll . 

Here, we investigate how the relative strengths of various heat sources affect the temperature 
distribution in the resonator. By performing numerical modeling, and measuring the locked 
laser's frequency change as a function of various thermal parameters, we gain an understand- 
ing of the mechanisms limiting performance and the steps necessary for their mitigation. In 
particular, we find that the dual-mode technique leads to unstable temperature gradients in the 
resonator that pull the frequency through thermal strain. Despite these limitations, the tech- 
nique reduces the effective temperature coefficient of the WGMR reference cavity by more 
than an order-of-mag nitude-from 8.88xl(T 6 K" 1 to 1.33xl(T 7 K _1 . We achieve an optical 
frequency instability of 1.29x 10 at one second, and a long-term drift of 21 kHz/hr. 

2 Dual-mode stabilization 

As described in detail previously [9 1, the thermorefractive coefficients of orthogonally polarized 
modes in a birefringent resonator are significantly different. Therefore, the frequency spacing 
between ordinary and extraordinary modes depends heavily upon temperature. The dependence 
is given by d(Af)/dT = — /()(<% — oci^), where 0£^ are the thermorefractive coefficients 
of the ordinary and extraordinary polarizations, and fo is the optical carrier frequency. For 
magnesium fluoride (MgF2) at our experimental temperature of 40 °C and fo = c/1560 nm, 
the differential tempco is 67 MHz/K. With a multi-mode resonator such as this, one can easily 
identify an appropriate pair of modes with spacing in the tens-of-megahertz range because, 
for each polarization, there are several tens of high- Q modes within the free spectral range of 
10.6 GHz. With the narrow linewidths in a MgF2 resonator, a frequency spacing measurement 
on the Hz level is readily achieved, yielding nK sensitivity. 

2.1 Experimental setup 

Our implementation of the dual-mode stabilization scheme is shown in Fig. Q] and described 
as follows. A narrow linewidth fiber laser is coupled via an angle-polished fiber lfl3l into 
3.25 mm radius, z-cut resonator made from excimer-grade MgF2. Mode linewidths are in the 
range of 40-80 kHz, corresponding to quality factors of 2.4-4.8 x 10 9 , depending on the mode 
order, polarization, and loading conditions. Laser light passes first through an acousto-optic 
modulator (AOM), and then an electro-optic phase modulator (EOM), before being coupled into 
the cavity at a polarization angle that can excite both modes. Two modulation frequencies are 
applied at the EOM: a 1 MHz signal for generation of Pound-Drever-Hall error signals 1141 . and 
a second signal at the frequency spacing Af between a chosen pair of orthogonal modes. Thus, 
the laser is frequency locked to the ordinary mode, while the modulation sideband coincides 




Figure 1: The experimental setup, (a) Schematic of the optical system and stabilization loops. 
Each proportional-integrator (PI) loop has both a slow and fast actuation branch. Polarization 
controller (PC) allows for control of excitation ratio between the two modes, (b) Photo of the 
resonator-TEM assembly, (c) Scope traces of the detector and PDH error signals as the laser is 
swept over the orthogonal mode pair; modulation frequency is 1 MHz. 

with the extraordinary mode, providing an error signal for temperature stabilization. For these 
experiments, the error signal has a noise-limited temperature discriminant of 140 nK Hz -1 / 2 . 
Thermal feedback (unity gain bandwidth w 1 Hz) to the resonator is provided by a 4 x 4 mm 2 
thermoelectric module (TEM) glued directly to the disc; an optional fast thermal loop actuates 
on the in-coupled laser intensity via the AOM. The resonator-TEM assembly is mounted in a 
small (2.5 x 2.5 x 2.5 cm 3 ) brass oven, whose temperature is measured with a thermistor and 
stabilized by feeding back to a resistive heater. 

The stability of the resonator coupling is critical to the frequency stabilization and its prac- 
tical use, so we paid special attention to the fiber coupler mounting scheme. The small oven 
in which the resonator is mounted has a base that is machined from a single block of brass. It 
has a flexure arm where the fiber coupler is attached. Initial alignment of the fiber coupler is 
performed with a 3D piezo stage, after which the fiber is glued to the arm. A micro screw is 
used on the flexure arm for subsequent small gap adjustments. This mounting scheme provides 
stable operation for periods as long as months without any need for adjustment and with no 
noticeable change in the coupling condition. 

The oven-resonator-coupler setup is placed in a cylindrical vacuum chamber that is evacuated 
to the milli-torr level. Temperature of the vacuum canister itself is not actively controlled. Use 
of the vacuum chamber eliminates convective heat transfer and provides isolation from changes 
in atmospheric pressure and humidity. 

We measure and characterize the system by examining the in-loop error signals and by mea- 
suring the beatnote between the locked laser and a self-referenced frequency comb. The comb 
is stabilized to a hydrogen maser that has a fractional instability less than 1 x 10~ 13 t~ 1//2 from 
1 to 10,000 s. 

2.2 Simulation 

The dual-mode stabilized resonator system has complex thermal properties and dynamics. To 
understand the various factors influencing mode volume temperature, as well as the overall 
thermal distributions in the resonator, finite element method (FEM) modeling proves useful. 
Two-dimensional, axially symmetric modeling of the resonator-TEM assembly is performed 



using a commercial FEM package (COMSOL). Material properties are taken from the manu- 
facturer datasheets or material property databases; model geometries are based as closely as 
possible on the physical system. Screenshots of the FEM model are shown in Fig. [2] 

The objective is to determine the frequency stability of a WGMR whose mode volume tem- 
perature is stabilized to the nK level via the dual-mode technique. Our model accounts for the 
following heat transfer mechanisms: radiative coupling between resonator and oven, optical 
heating of the mode volume, heat pumped into/out of the resonator by the thermoelectric ef- 
fect, and thermal conduction through the TEM. (We assume convective transfer is negligible in 
the vacuum chamber.) Here, we are primarily concerned with dynamics within the TEM feed- 
back bandwidth, i.e., those with times-scales >1 s. We calculate the thermally-induced strains 
on the resonator's radius r and how they pull the frequency through A///o = —Ar/r. The 
model does not account for thermorefraction because the mode volume temperature is taken to 
be stable, as indicated by the residual error signal of the temperature stabilization loop. 

Radiative heat transfer between one body completely enveloped by another can be estimated 

as na 

^- 1/e , + ,,M 2 (ife- 1 ) ' r '-^- <» 

where e is the emissivity, A the area, T the temperature, the subscripts indicate inner (1) and 
outer (2) bodies, and a is the Stefan-Boltzmann constant. For various values of temperature 
difference between resonator and enclosure, q mi j is calculated using Eq. (Q]i and then applied as 
a boundary condition on the resonator surface. 

Laser light circulating in the mode is converted to heat through material absorption. For a 
laser locked on resonance (zero detuning), the heat generated by this absorption can be esti- 
mated through |fl6l 

qiaser=lT]-p: (2) 

Slabs 

where / is the input laser power, r\ is the coupling efficiency, and Q (Q a b s ) are the observed 
(absorption limited) quality factors. This heat is included in the model as a volume heat source 
(units of J/m 3 ) applied to the mode volume (region shown in pink in Fig. |2a)). The mode 
volume cross-sectional shape is an ellipse with major (minor) axes of 20 (10) fim, which ap- 
proximates the distribution of the fundamental mode as calculated by a separate, full vectorial 
FEM simulation as described in ifTTll . 

Heat generated in the TEM (qrEM) is modeled as a volume source distributed over the top 
alumina plate. The bottom of the TEM is held at the enclosure temperature. Thus, we can model 
how the three heat sources q ra d, qiaser, and qrEM contribute to the mode volume temperature. 
An initial steady-state solution is found assuming reasonable values of in-coupled laser power, 
with qrEM = Qrad = 0, yielding values of disc radius and mode volume temperature. The 
dual mode condition is simulated by demanding that subsequent changes in q rai i or qi ase r are 
compensated for by a change in qrEM, maintaining the mode volume temperature at its initial 
value. The resulting strains (changes in radius) are converted to frequency, allowing us to model 
how a change in enclosure temperature or laser intensity that occurs while the dual-mode loops 
are closed will affect the frequency of the resonator. 

3 Results 

The dual-mode scheme measures and stabilizes the temperature in the region occupied by the 
optical modes. Our loop guarantees that any uncontrolled (out-of-loop) change in thermal 
transport parameters will result in a compensating change in the level of thermal feedback 
actuation. As the relative strengths of these thermal parameters changes, the system moves 



along lines of constant mode volume temperature. These isotherms do not, however, represent 
lines of constant resonance frequency because the central portion of the resonator changes 
temperature and expands/contracts. The isotherm slope, i.e., the ratio of change in optical 
frequency to the change in a system parameter, represents an effective coefficient of frequency. 



3.1 Environmental temperature dependence 

We simulated the dual-mode resonator's effective tempco by solving the FEM model for various 
temperature differences between resonator and enclosure (T\ and T2 in Eq. ([TJ). The condition 
T\ = T2 means the average resonator temperature is in equilibrium with the oven and represents 
the instant at which the dual-mode loops are closed. This q ra[ { = situation is solved first. 
Next, a small drift in T2 is imposed, q ra j is calculated, applied as a boundary condition on the 
resonator, and the necessary value of qjEM needed to return the mode-volume temperature to 
its initial value is found - this represents the feedback action of the dual-mode loop. 
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Figure 2: FEM simulation, (a) Model setup showing the mesh, materials, a detail of the mode 
volume (pink region), and the line along which temperature and strain are plotted (dashed or- 
ange). BiTe: bismuth telluride. (b) Temperature map results showing how different gradients 
exist for different relative strengths of heat sources. Plots of temperature (c) and radial de- 
formation (d) along the disc radius for T2 = T\ (solid red) and T2 = T\ + 50 mK (dashed blue), 
showing how the latter distribution is cooler in the center, contracting the disc radius by «6 pm. 

The result is illustrated in Fig.|2b), where the different thermal distribution resulting from a 
change in enclosure temperature is clearly visible. At T2 = T\, no radiative heating of the disc 
occurs; the gradient forms as optically-generated heat is conducted out through the TEM. If 
the oven temperature drifts up by 50 mK, radiation additionally heats the top and outer portion 
of the disc. Maintaining the mode volume temperature - out to the nK digit in these simula- 
tions - requires that the TEM cool the central-bottom portion of the resonator, contracting it 
and increasing the common-mode optical frequency. The plots in Fig. [2] show the calculated 
distribution of temperature (c) and deformation (d) along the dashed orange radial line in (a). 
The temperature distribution is normalized to the mode volume temperature (Ty m ), and the dif- 
ference in strains at the endpoint (R = 3.25 mm) yields the frequency shift. The calculated shift 
was found to be linear across five values of oven temperature drift spanning half a degree. A 
linear fit to these simulated data yields an effective tempco of 6.9 kHz/mK, corresponding to a 



fractional coefficient of 3.6 x 10~ 8 K _I . We note that these values are comparable to the linear 
thermal expansion coefficients of common ULE glasses. 

The effective tempco was also measured experimentally. With the laser locked and the dual- 
mode stabilization loops engaged, we increased the setpoint on the oven's temperature con- 
troller in 12 mK steps every 360 seconds (eventually returning it to its initial value). The TEM 
current (converted to heat via values on the TEM's datasheet), locked laser's optical frequency, 
and in-loop temperature error signal were simultaneously recorded during the run, and are pre- 
sented in Fig. [3] In agreement with the simulation, increased radiative heating of the disc results 
in compensatory cooling from the TEM, leaving the mode volume temperature the same but the 
optical frequency higher. A linear fit of the optical frequency as a function of change in enclo- 
sure temperature gives an effective tempco, with the experimental value being 25.5 kHz/mK. 
The magnitude of this value is 67 times less than the thermal-expansion determined tempco of 
-1.71 MHz/mK, and it corresponds to a fractional coefficient of 1.33 x 10 K . 




Figure 3: Effect of environmental temperature: Time series (a) showing shift in optical fre- 
quency (Af opt ) and heat delivered to the resonator by the TEM as enclosure temperature 
changes in 12 mK steps. Inset shows mode volume temperature from the in-loop error signal; 
aside from spikes at the temperature steps, it remains extremely stable. The multi-dimensional 
relationship between parameters is shown in (b). T2 — T\ is the drift in enclosure temperature. 
Triangles are data from (a) averaged over 360 s bins. The isotherm is a fit to these data (solid 
blue line), while the slope of its projection onto the temperature-frequency axis (dashed red 
line) yields the effective frequency coefficient of 25.5 kHz/mK. 



3.2 Circulating power dependence 

Despite MgF2's extremely high transmission at 1560 nm, absorbed optical power is constantly 
converted to heat in the mode volume. Thus, changes in circulating optical power will desta- 
bilize the optical frequency at any time scale. Here, we investigate the effect of slow (within 
the TEM loop bandwidth) circulating power changes. These slower changes may be harder to 
control since they result from drifts in polarization or coupling strength, while faster fluctua- 
tions are in principle corrected for by the intensity branch of the dual-mode loop. Similar to the 
case of environmental temperature, we find a multi-parameter dynamic whereby a change in 
circulating power results in a change in TEM correction signal, thermal gradient, and resonator 
frequency. 

Using the same FEM model described above, we calculated the frequency shift as a function 
of change in circulating optical power. For these simulations q ra( i was kept at zero. Heat 
was allowed to leave the resonator radiatively through its surface (to a boundary at infinity) 
and via conduction through the TEM. For a given change in circulating power, the resulting 



change in qjEM necessary to maintain the mode volume, change in thermal gradient, and shift 
in frequency was calculated. Using Eq. (J2j and reasonable values of / = 100 }X W at the coupler, 
r; = 0.07, and that the disc's quality factor is absorption limited, the simulations suggest that a 
5% reduction in laser power will yield a frequency shift of -103 kHz. 

An experimental measurement of this isotherm slope was also made. With the laser inten- 
sity branch of the thermal loop open, we changed the optical power sent to the resonator by 
varying the amount of light transmitted through the AOM. In agreement with the simulations, a 
decrease in circulating power causes an increase in TEM heating to maintain the mode volume 
temperature. The new thermal gradient is warmer in the disc center, increasing r and decreasing 
the optical frequency. The experimentally measured frequency pulling was -15 kHz for a 5% 
reduction in intensity, corresponding to an effective coefficient of 43 kHz//iW of circulating 
optical power. Based on this result, a stability of ~47 pW would be needed to reach a bench- 
mark fractional frequency instability of 10~ 14 . This requirement on laser intensity stabilization 
is a readily achieved 7 ppm, however we stress that the requirement is on circulating power, so 
coupling strength and polarization angle (which changes the power ratio between modes) also 
require ppm-level stability. 

The discrepancy between simulated and experimentally determined temperature coefficients 
results from the heavy dependence on setup details such as boundary conditions at the domain 
interfaces and values of material properties; assumptions made in using Eqs. (1) and (2) could 
also lead to errors. Experimentally, it is difficult to isolate a single variable as can be done with 
simulation, so various other destabilizing factors can play a role. However, the simulations 
and experiments agree very well qualitatively, provide meaningful insight, and with refinement 
could agree quantitatively as well. 

3.3 Dual-mode frequency stability 

To evaluate the performance gained from the dual-mode temperature stabilization technique, 
we made optical comparisons to a frequency comb. After the AOM, a portion of the stabilized 
laser light was picked off and beat on a photodiode with a bandpass-filtered span of the comb. 
For several thousand seconds, the beatnote was counted with 10 ms gate times. Results of two 
typical data runs are shown in Fig. |4] The lowest instability of 1.29 xlO was reached at 
700 ms, while a linear fit to the 6000 second time-series frequency excursions has a slope of 
21 kHz/hr. 

In Fig. |4] we show data for both the dual-mode stabilized and "free-running" resonator, 
where the latter configuration is identical except that the temperature feedback loop is open. 
Stabilities below 0.1 s are identical since the dual-mode locking loop has little or no gain at the 
corresponding frequencies. The bump just below 0.1 s has unknown origin but is likely due 
to laser noises. Above 0.1 s the free running system drifts linearly with time. The dual-mode 
system reaches the stability floor at 1 s and then starts a random walk with time largely due to 
temperature gradient variation as discussed in section 3.1. 

Because of these gradient effects, the effective tempco is not completely eliminated. Overall 
frequency stability of our dual-mode setup therefore remains limited by temperature fluctua- 
tions of the brass oven in which the resonator is mounted. An analog servo controls the oven's 
temperature, and the loop oscillates with a period of «2000 s. The amplitude of the oven's 
temperature oscillations is 2-3 mK, consistent with the frequency excursions seen in the inset 
of Fig.|4]for both the free running and dual-mode tempcos of -1.71 MHz/mK and 25.5 kHz/mK, 
respectively. However, Fig.|4]demonstrates that for a given environmental instability, the dual- 
mode loop yields improvements of more than an order of magnitude within the dual-mode 
bandwidth, i.e. for T > 1 s. 

Frequency drifts caused by changes in laser power or coupling conditions remain below the 
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Figure 4: Overlapping Allan deviation of a laser locked to the WGMR cavity, as compared 
to the frequency comb, under nominally identical conditions except that the dual-mode loop 
is closed (blue triangles) or open (red circles). Significant improvement is observed for time 
scales > 100 ms. Inset: Time sequences of the data run showing frequency excursions for 
the dual-mode (blue, left scale) and free-running (red, right scale) resonator, showing periodic 
oscillation from the enclosure's analog temperature controller. 

limit set by environmental instabilities in our setup. Therefore, better temperature stabilization 
of the oven would yield better overall performance. To achieve an even higher frequency stabil- 
ity, temperature gradient effects must be reduced. Nevertheless, the experimentally determined 
effective dual-mode tempco is less than an order of magnitude larger than that of ULE cavities, 
and although our efforts to stabilize the environmental temperature are modest, the long-term 
drift rate of 21 kHz/hr is among the best stability results achieved with a whispering gallery 
mode reference cavity. 

4 Conclusion 

We have investigated the frequency stability of a laser locked to a dual-mode WGMR refer- 
ence cavity. Because the temperature distribution in the resonator depends on several thermal 
parameters simultaneously, we find that a precisely stabilized mode volume temperature does 
not guarantee corresponding frequency stability. Nonetheless, we calculate that the dual mode 
technique can reduce the temperature coefficient to 3.6xl0~ 8 K _1 , comparable to the linear 
coefficient of ULE. In addition to environmental temperature, the level of circulating optical 
power also plays a significant roll. We experimentally demonstrate a reduction of temperature 
coefficient to 1.33 x 10~ 7 K _1 , a frequency instability at one second of 1.29x 10~ 12 , and a drift 
of 21 kHz/hr over long time periods. These results represent progress towards using WGMRs 
as optical references for frequency metrology applications by partially suppressing the high 
temperature dependence of WGMRs, thus relaxing the requirement on environmental stability. 
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